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ABSTRACT 



Time histories of the development of the magnetosphere boundary 
following impact of a sharp) -fronted plasma cloud and of the compression 
of the magnetosphere following inpact by a discontinuity in the solar 
wind are calculated. Results are presented for initially planar discon- 
tinuities having their normals either aligned with or inclined 30 ° to the 
wind direction. Although several minutes are generally required for a 
discontinuity in the solar wind to sweep past the geophysically significant 
portion of the magnetosphere, it is found that any individual element of 
the magnetosphere boundary usually conpletes the major part of its adjust- 
ment to the new conditions within a period of about a minute. 
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1. INTRODUCTION AND FORMULATION OF PROBLEM 

Although the basic concepts fundamental to understanding how the geo- 
magnetic field carves a hollow in a rapidly streaming solar plasma derive 
principally from a long series of pioneering studies of transient phenomena 
by Chapman and Ferraro (see CHAPMAN, 1963 , for a recent review), most recent 
studies have been concerned with the corresponding steady-state problem 
(see BEARD, 1964, for a recent review). Since the latter problem is now 
becoming relatively well solved and the results generally agree with 
measurements made in space (see particularly NESS, SCEARCE, and SEEK, 1964), 
we return to two unsteady -state problems similar to those originally con- 
sidered by Chapman and Ferraro over 30 years ago. We consider first that a 
huge isolated cloud of plasma advances toward and engulfs the earth, and, 
second, that a discontinuity in momentum flux density in an otherwise steady 
solar wind sweeps past the earth. The latter more realistic possibility 
may be associated either with a shock wave across which both the density 
and velocity are discontinuous, or with a contact surface acrdss which the 
density, but not the velocity, changes. 

t 

We shall consider only the geomagnetic equatorial cross ‘ section of the 
hollow, and that the velocity vector V of the undisturbed incident stream 
lies in this plane. We introduce rectangular Cartesian coordinates with 
the origin at the center of the earth, the xy plane lying in the geomag - 

i 

netic equatorial plane, and the y ( axis directed upstream as illustrated 
in Fig. 1. I 

The Magnetic field at the magnetosphere boundary will be denoted by 
B s , and the angle between y and the normal n to the boundary by i|r. 
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The dynamical balance of magnetic and fluid pressures on the two sides of 
the moving boundary provides the following boundary condition at the surface 
of the hollow (FERRARO, 1952, i960) 

B s 2 /8jt = Kmn[ ( V - v) • n] 2 = Kmn[ V cos ^ - v n ) (l) 

where m and n are the mass and number density of ions in the incident 
stream, v is the velocity of a segment of the boundary, and 
v n = v • ft = v cos + is its component normal to the boundary. Cosine ^ 
can moreover be expressed as follows : 


cos 


= dx/(dx 2 + dy 2 ) 


1/2 


( 2 ) 


where time is to be held constant in the evaluation of the partial derivatives. 
The coefficient K depends on the nature of the interaction between the 
particles of the solar wind and the boundary. Specular reflection, for which 
the appropriate value is K = 2, has been preferred historically, but recent 
developments (see, e.g., SPREITER and JONES, 1963) suggest that a hydromag - 
netic fluid model, for which K = 1, is to be favored. 

Following FERRARO (1952), the unknown magnetic field B s at the 
magnetosphere boundary will be approximated by twice the tangential compo- 
nent of the earth's magnetic dipole field. Thus, 

B s = 2B 0 (a/r) 3 (3) 

where B 0 = 0.312 gauss is the equatorial value of the earth's surface mag- 

„ = I 1 => 9,1/2 

netic field, a = 6.37x10 s cm is the radius of the earth, and r = (x* + y*) 
is the radial distance from the center of the earth. 
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Solving equations (l) and (3) for v n /V and selecting the proper sign 
so that 


in the final 


cos 

steady state when 


v n 

— = cos 


+ = Bs/SrtPst 

v n = 0 yields 

♦ - 52 

(2*Pst> l/2 
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(5) 


where p g ^ = KxnnV 2 represents the stagnation pressure at the apex of the 
hollow. Cartesian components of the normal velocity of a surface element 
are thus 

- v n sin ♦ , ^ = -v n cob * (6) 


Given a curve or set of coordinates for the initial equatorial trace of the 
geomagnetic boundary, it is possible to integrate equation (6) numerically 
and determine the corresponding boundary coordinates for a slightly later 
time. Successive application of this process enables the calculation of 
the boundary coordinates for all later times. 

t 

The above equations may be expressed alternatively in dimensionless 
form by introduction of the distance p Q from the center of the earth to 
the apex of the steady-state hollow (v n = 0, t = 0) as a unit of length, 
and the time t 0 required to travel a distance p Q at speed V as a 
unit of time; that is, i 


a(B 0 2 /2np, 


st 


> 1/s . 


- °o/ T 
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Thus, equations (5 ) and (6) can be rewritten as follows: 
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where 

i = x/p 0 > H - y/p 0 > P = r / P Q > t = t/r 0 


(9) 


2. INTERACTION WITH AN ISOLATED PLASMA CLOUD 

Over 30 years ago, CHATMAN and FERRARO (1931) described how an approxi- 
mately parabolic hollow would be formed in an isolated and initially flat- 
faced plasma cloud as it approaches and passes beyond the earth. Although 
their sketches have been reproduced widely, and FERRARO (1952, i960) has 
more recently integrated equation (6) with ^ = 0 (see also CHAPMAN, i960, 
for a review), the time history of the motion of the remainder of the front 
of the plasma cloud has never been computed. We shall discuss this problem 
first because of its relative simplicity, conceptual value, and historical 
interest. 

We consider that a uniform plasma cloud of density mn sweeps - toward 
the earth with velocity V and that its leading surface is planar at time 
t = 0 when at a normal distance r^ from the center of the earth. The 
distance r i is selected so as to be sufficiently large that the previous 
deformation of the boundary as it approaches r^ from infinity is small, 
but this is actually immaterial to j the calculations since we are only 

concerned with the behavior of the; plasma boundary at times subsequent to 

1 

I 

the initial time. The result may be presented either in convenient dimen- 
sional form or more universal dimehs ionless form. We have elected the first 
alternative, since the scaling laws are very simple, and have selected the 
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following representative values for the pertinent variables : 

V = 500 km/sec = 5xl0 7 cm/sec; n = 5 protons/cm 3 ; m = 1.67xl0“ 24 gmj and 
K = 1. Results are shown in Fig. 2 for the case in which the initially 
planar front face of the plasma cloud is parallel to the x axis at a 
distance r^ = 20 earth radii. The corresponding results are shown in 
Fig. 3 for a case in which the initial planar face is inclined 30° to 
the x axis. The dashed line shows the final steady -state configuration 
of the boundary defined by the integral curve of equation (4) that passes 
through p = 1 at x = 0. The latter solution was indicated first in 
parametric form by FERRARO (1952), and evaluated explicitly in numerical 
form by BEARD (i960) and SPREITER and BRIGGS (1961). 

The results in Figs. 2 and 3 are free of surprises and are quite con- 
sistent with the early qualitative descriptions of Chapman and Ferraro. 
Quantitative features of particular interest are the rapidity with which 
the boundary approaches its steady-state form as the plasma front sweeps 
by, the smallness of effects of obliquity of the plasma front, and the 
tendency to form waves that drift downstream along the boundary as it nears 
the steady-state cohfiguration. The latter tendency is consistent with 
analytic results that can be derived for small scale waves on the magneto- 
sphere boundary. 

The results described above can be converted by a simple relabeling 
to those for other values for the pertinent variables by use of the following 
relations derived from equation 

■ L 2 /i i = x 2/ x i 3 = ^stj/^l 
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(9) 
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where subscript 1 refers to the values used to calculate the results shown 
in Figs. 2 and 3> and subscript 2 refers to any other set of values. 

3. INTERACTION WITH A DISCONTINUITY IN MOMENTUM FLUX DENSITY IN THE SOLAR WIND 


We consider next that a steady solar wind has existed for a sufficiently 
long time (of the order of 10 minutes on the basis of the preceding results ) 
that the boundary between the geomagnetic field and the solar plasma is well 
established in the steady -state configuration associated with a value p g ^. 

61 

for the stagnation pressure. We consider further that a discontinuity in 
momentum flux advances through or with the solar wind plasma with velocity 
V d , and that the discontinuity surface is planar over lateral distances 
larger than the magnetosphere before it begins to interact with the magnetic 
field of the earth. Behind the discontinuity surface, conditions are consid- 
ered to be steady and uniform, but are characterized by a larger stagnation 
pressure p .. Ultimately, therefore, the magnetosphere boundary will assume 
a configuration geometrically similar to that possessed originally, but 
smaller in the ratio (p 8 ^. /p s ^) . We seek here to calculate the' time 

history of the collapse of the magnetosphere from the initial to the final 
configuration. 

The governing differential equations are the same as in the preceding 
case, but the initial conditions are different. Thus, the initial shape of 
the magnetosphere boundary is given by the solution of equation (4) with 


p st = P s t * A-fc an interval of time At after the discontinuity surface 
passes the magnetosphere nose, thi most upstream portion of the magneto- 
sphere extending over a distance V^At in the y direction will have 
been exposed to the greater pressure P 8 ^ behind the discontinuity 
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surface and set into motion toward its final equilibrium position. The 
remainder of the magnetosphere boundary has not yet experienced "the enhanced 
pressure, however, and remains fixed in its initial configuration. A gen- 
eral situation 'that may exist at any time after the discontinuity surface 
has passed the initial position of the magnetosphere nose is illustrated 
in Fig. 4. The shape of the magnetosphere boundary and the discontinuity 
surface at a s m all time interval later is calculated as follows. First, 
the part of the discontinuity surface exterior to the initial magnetosphere 
shape simply advances with velocity and remains straight and parallel 

to its initial direction. The remainder of the discontinuity surface coin- 
cides at any moment with a portion of the magnetosphere boundary, and 
advances with velocity components given by equation (6) using the values 
and P s tb for V and p s f The displacement of this portion of the mag- 
netosphere boundary during a small time interval can then be determined 
numerically in the same manner as in the preceding case. Since each bound- 
ary element moves normal to itself in this calculation, there always remains 
small unconnected gaps separating the ends of the moving portion of the mag- 
netosphere boundary and the straight portion of the discontinuity surface 
exterior to the magnetosphere. The necessary connections have been made 
with a parabola for the first time interval after impact of the discontinuity 
surface and the magnetosphere nose, and with straight lines for all subsequent 
times. The remainder of the magnetosphere boundary remains, of course, in 
its initial configuration until the arrival of the discontinuity surface. 

Results are presented in Figs. 5 and 6 for two specific cases in which 
the initial planar discontinuity surface is inclined at 0° and 30° to the 
x axis. They are for a discontinuity surface which advances with a velocity 
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V d = 5xl0 7 cm/sec and is characterized by a fourfold increase in p s ^. from 
p s t a = 2.09xl0 " 18 dynes/cm 2 to p st ^=hp s -t • This value for P B t a could 
represent, for instance, a flux of 5 protons/cm 3 traveling at 500 km/sec 
with K = 1, or. any other combination of values for which the product 
p s -(. = (KmnV 2 ) a is the same. The initial and final steady -state solutions 

are indicated by dashed lines. As in the previous examples, the boundary 
is found to collapse quite rapidly to its final configuration. Although 
several minutes are required for the discontinuity surface to sweep past 
the forward part of the magnetosphere of greatest geophysical interest, 
each element of the boundary essentially completes its movement from 
initial to final position in less than a minute. 

The following relations can be used to convert the results in Figs. 5 
and 6 to those associated with other values of the stated parameters; 



where the notation with respect to subscripts 1 and 2 is the same as in the 
case of equation (10). It is required that the ratio P s t^/Pst a = ^ as 
stated above, however. Results for other values for this ratio could be 
computed in the same manner, but there is little reason to believe the 
results would differ substantially from those presented here. 
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